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Abstract 

A close inspection on the 3D hydrogen atom Hamiltonian revealed 
formal eigenvectors often discarded in the literature. Although not 
in its domain, such eigenvectors belong to the Hilbert space, and so 
their time evolution is well defined. They are then related to the ID 
and 2D hydrogen atoms and it is numerically found that they have 
continuous components, so that ionization can take place. 
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1 Introduction 



In order to clearly state the question addressed here, it is important to re- 
call some points of the mathematical foundation of observables in quantum 
mechanics. There will be two main contributions, one related to some weak 
solutions of the Schrodinger equation and other to dimensional interpreta- 
tions. 

The problem of finding the correct self-adjoint extension describing the 
quantum (Schrodinger) operator corresponding to a physical model can be 
subtle and difficult. Usually the physicist has a clear expression for the 
operator, an unbounded one acting in a Hilbert space 7i, but it is not obvious 
which domain should be taken (some general references for what follows are 

BEIE!). 

Let ('0, 0) denote the inner product in if T is a linear operator acting 
on its dense domain domT C 7i, then to represent a physical observable it 
is necessary that T is hermitian, i.e., 

{Tip,(f)) = {ip,T(f)), V'?/', G domT. 

However this condition is not enough to guarantee that T has real spectrum 
and the time evolution it generates is unitary; the right condition is self- 
adjointness. The domain of its adjoin T* is 

domT* = {^eH:37]en with (r/, 0) = (^,T0), V0 G domT}, 

and for ^ G domT*, one has T*^ = t]. It follows that T* is well defined if 
domT is dense in 7i, and T is hermitian if, and only if, T* is an extension of T. 
The operator T is self-adjoint if T = T*. Notice also that (often) for bounded 
operators the distinction between hermitian and self-adjoint operators does 
not exist. 

As already mentioned, usually T is hermitian with dense domain, and one 
asks if it is also self-adjoint or has any self-adjoint extension; such extensions 
are the candidates for the operator describing the related physical observ- 
able. A nice situation that often occurs, in particular for the Hamiltonian of 
the Hydrogen atom (and other atomic systems as well), is that T is essen- 
tially self-adjoint, i.e., it has just one self-adjoint extension and the physical 
operator is well determined. However, there are situations where there are 
infinitely many self-adjoint extensions and each one should correspond to a 
different physical circumstance; the choice is a physical one, not on math- 
ematical bases. Even worse, some hermitian operators have no self-adjoint 
extensions! 

The standard example of such framework is the momentum operator P = 
—i-^ for a particle in a box [0, 1]. In this case Ti. = T^[0, 1], it is natural to 
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take domP as smooth functions ip (zTi such that ^/'(O) = = (so that 
the particle remains confined to the box); the self-adjoint extensions of this 
hermitian operator are Pq,, where a is a complex number with |a| = 1, and 
all elements of domPo, satisfy il^iX) = aip{0). 

It is worth remarking that if T is hermitian and domP = Ti, then T is 
bounded, so that in general such domain questions are not avoidable. These 
interesting problems are well explored in the literature, and as additional 
references see ji] and for applications to the one-dimensional hydrogen atom 
see 0IE]. 

Nevertheless, there are some delicate issues in the mathematical founda- 
tions of quantum mechanics that seem not yet exploited from the physical 
point of view. The main goal of this work is to discuss one of such issues and 
relate it to a physical situation. 

Recall that a self-adjoint Hamiltonian operator H generates a time evo- 
lution ipit) = U{t,0)ip = e-'^^ip, which is a solution of the Schrodinger 
equation 

Since U{t, 0) is a family of unitary operators, for any time t its domain is 
the whole Hilbert space 7i, so that it is meaningful to consider [/(t, O)^? for 
If & 7i but with (p ^ domH, i.e., the time evolution is not restricted to the 
domain of H. Sometimes U{t,0)if, for (f ^ domH, is called a weak solution 
of the Schrodinger equation. 

An unusual situation will be presented. A self-adjoint operator H with 
dense domH izTi = L^(IR^) will be considered, vectors S G ?i not belonging 
to its domain will be given, although they are pseudo-eigenvectors of if, that 
is, 

H^ = AhS, (1) 

for Ah G 1R. Some numerical calculations will indicate that S has a nonzero 
component in the continuous subspace of H, so that the naive time evolution 
built from gives an incorrect answer. It will be argued that such solu- 
tions are related to the same model in smaller dimensions. Furthermore, the 
physical system in question is one of the most celebrated models in quantum 
mechanics, the three dimensional (3D) hydrogen atom. 



2 Pseudo-Eigenvectors as Weak Solutions 

The hermitian Hamiltonian of the 3D hydrogen atom is 

Ho = -—A , domHo = C^i'R^) C L^JR^), 
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where /i is the electron mass, e its electric charge and C^(IR^) denotes the 
set of smooth functions with compact support. This operator is essentially 
self-adjoint and its unique self-adjoint extension Hh-, the 3D hydrogen atom 
Hamiltonian, reads 



H 



H 



A 

2/i r 



dom Hh = H^(JR^ 



(2) 



with if^(lR^) denoting an appropriate Sobolev space; in particular, if^(lR^) 
is a subspace of L^(IR^), it is also the natural domain of the free particle 
Hamiltonian and all its elements are continuous functions ^2|. 

The usual spectral analysis of Hh can be performed and its well-known 
eigenvalues 



/ie 



n > 1, 



2/1 

can be found. Recall that the closed subspace Hp generated by its eigenvec- 
tors is named the point subspace of Hh and its orthogonal complement Hac 
is a nontrivial subspace (i.e., it has nonzero elements) and named the abso- 
lutely continuous (or scattering) subspace of Hh- Physically, the members 
of TCp are the bound states while the elements of Hac describe the ionizing 
atomic states (this interpretation follows, for instance, by the RAGE Theo- 
rem laiH]). 

The eigenvalue equation for the 3D hydrogen atom Hamiltonian is sepa- 
rable in spherical coordinates r>0, O<0<7r, O<0< 27r, and by taking 
the standard representation 



i?(r)e(^)$(0). 



the equation for 0(6') is given by [TT] 



m 



sin^ 9 



+ 1) 6 = 0, 



(3) 



(4) 



with m and i > being integer constants. For each i value one has —i < 
m < i. 

Consider first the particular case £ = 0; it follows that m = and 
reduces to 

1 d ( dQ\ 

The usual normalized solution of this equation is 6o,o(^) = 0z=o,m=o(6') = 
1 / However, there is also the additional solution (that will play a major 
role here) 



eo,o(^) = — In 

TT 



tan 



(6) 
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This is just one instance of additional solutions S,i^m of for £, m as above; 
such solutions are Legendre function of the second kind IHl HDl- The ^^^^ 
solutions have been discarded in the mathematical literature since they are 
not continuous at 6^ = and 9 = n, and so via (jHI) they do not generate 
elements in the domain of Hh', and discarded in the physical literature ^] 
by arguing they are not bounded functions. 

By taking the usual radial Rn/{r) and azimuthal $m(0) = solu- 
tions for the 3D hydrogen atom, set {£ < n) 



2 \3 

nao/ 47m [(n + £)!]■'_ 



exp 



(-) ^» 



2£+l 



nao/ \nao/ \naQ 
so that one gets the standard eigenfunctions (here restricted to m = 0) 

^„,,,o(r, e, 0) := Fn,i,o{r, 9, 0)0^,o(^), 



and now the additional ones 



S„,£,o(r, 0, (p) ■= Fn/fl{r, 9, 0)6,0(6'), 

where oq = h^/{^e^) denotes the Bohr radius and L^i^l are the Laguerre 
polynomials. 

For m 7^ the probability density generated by diverges (recall the 
Jacobian is r^sin6'), i.e., 

r sme\i,,^^,{e)\^de = ^■ 
thus the corresponding functions 'E.n,i,m^o{r, 9, 0) do not belong to the Hilbert 
space L2(]R3). Hence, it is meaningless to talk about their time evolution 
even as weak solutions of the 3D hydrogen atom Schrodinger equation. 

However, for m = the probability density generated by C,e,o does not 
diverge, i.e., by choosing appropriate constants it is found that 

Jo 

so that together with the R{r) and $(0) counterparts in Q, 6,0 generates 
elements S„,£,o of -L^(IR^), as given by the expression above. 

Notice that although En/fl does not belong to the domain of Hh, one 
formally finds 

Hn'^n/fl = ^n'^n/fl, = 2~2' < < "^J (7) 
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so that and A„ are pseudo-eigenvectors and pseudo-eigenvalues of Hh, 
respectively. Another point supporting the use of the adjective "pseudo" 
is that 'En,efl are not orthogonal to every \E'n,£',o (for instance, H„^o,o is not 
orthogonal to \E'n,^',o with odd £'). 

The question to be addressed now is about the time evolution of En/fi, 
which is well defined since 'En,e,o ^ L^(IR^) and so t/(t,0)S„_£o is a weak 
solution of the 3D hydrogen atom Schrodinger equation 

d 

ot 

Based on ((Zj) the naive expression for the solution 'En/fi{t) = e~'^^"^^^'^n,ifl is 

^nlfiit) = e ^^"^^/^^In/fi'i (8) 

such solution is not correct for o is not in the domain of Hh] if (|HI) holds 
then 'En/fiif) would be strongly differentiable and so one could conclude that 
Sn,£,o ^ dovciHH- Notice that (jHl) is correct if S is replaced by ^ . 

From the dynamic point of view it is important to check if q belongs 
to the point subspace Tip associated to Hh or if it has a component in the 
absolutely continuous subspace Tiac- In the latter case ionization can take 
place; in the former case these generalised eigenstates are written as super- 
positons of ordinary eigenstates, even if they do not belong to the domain of 
the operator, and so such solutions would be, in some sense, between linear 
combination of ordinary eigenvectors and the continuous space (but with no 
ionization) . 

In order to check if S„^£_o is generated by \E'n',^',m,, > 1, = 0, ■ ■ ■ , n' — 1, 
—I' < m < i', i.e., if 'En,e,o belongs to the point subspace of Hh, consider 

oo n' — l 
n'=l l'=0 

with Xne denoting the component of H„^ o in the continuous subspace Hac- 
Notice that clearly 'En/fi -L '^n',i',m if m 7^ 0. It was numerically found that 
IIXioll > O-S; as indicated in figure [T] for n = 1 (in figure [T] the values of 
P{N) are exact, since symbolic calculus was used); the parameter P(A^)^ = 

J2n'=i ^?=o^ l^^^il"^ is an approximation for 1 — ||xi,olP- Similar results were 
found for other values of i. 

Therefore one concludes that S„ £ have both nonzero point and continu- 
ous components, so that their time evolutions actually are not described by 
(jHI), but give nonzero probabilities HXn^lP of ionization (and then far from 
being bound states). 
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Figure 1: P as function of A^. 

3 Lower Dimensional Hydrogen Atom 

The fact that all q belong the Hilbert space raises the possibility of finding 
physical meanings for them; this section aims at discussing possible physical 
contents of these pseudo-eigenvectors. The first crucial remark is that for- 
mally H„^£ o has null azimuthal angular momentum (m = 0). The solution 
S„_o,o has also null total angular momentum (both i = = m), but there 
is a lack of rotational symmetry (see ©); this particular solution gives a 
clue on the physical interpretation. In fact, in comparision with ordinary 
eigenf unctions Sn/flir, 6) are elongated over the 2;-axis with a logarithmic 
divergence at 6* = and vr. Figure El shows the absolutely values of ^0,0(6') 
and Qofii^) as a function of 9, and Figure 01 a boundary surface of the 3D 
wavefunction Si 0,0, which is to be compared with ^^0,0,0 that has complete 
radial symmetry (its boundary surfaces are spheres centrered at the origin). 
Hence, there is a strong indication that S„^o,o are reminiscent of classical tra- 
jectories performing one-dimensional (ID) like motion, in agreement with its 
null angular momentum and lack of rotational symmetry. So, it is natural 
to relate such wavefunctions to the ID hydrogen atom, an interesting and 
controversial subject, popularized by the work of Loudon 12 published in 
1959. 

Loudon stated that the ID hydrogen atom was twofold degenerate, having 
even and odd eigenf unctions for each eigenvalue, except for the (even) ground 
state having infinite binding energy. Typically ID systems have no degener- 
ate eigenvalues, and Loudon justified the double degeneracy as a consequence 
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Figure 2: |.^o,o(^)| and |9o,o(^)| (dashed) as function of 6. 



of the singular atomic potential. Andrews jTH] questioned the existence of a 
ground state with infinite binding energy. Ten years later Haines and Roberts 
|14j revised Loudon's work and obtained that the even wave functions, with 
continuous eigenvalues, were complementary to odd functions, but such re- 
sults were criticized by Andrews fSl, who did not accepted the continuous 
eigenvalues. Gomes and Zimerman jTHj argued that the even states with fi- 
nite energy should be excluded. Spector and Lee [17] presented a relativistic 
treatment that removed the problem of infinite binding energy of the ground 
state. Several other works [IHlliniEniEIlESliniEHlEil (see also references 
therein) have discussed this (apparent) simple problem. 

The ID hydrogen atom has been used as a simplification of the 3D model 
in several theoretical and numerical studies 1211 12Z] • It is then interesting 
that Cole and Cohen |2H] and Wong et al. [211 have reported some experi- 
mental evidence for the ID hydrogen atom. The "quasi-lD" solutions H„ o,o 
are natural candidates to describe such experimental observations and may 
be relevant for an appropriate justification for the use of ID simplifications. 
Lastly, the ID eigenvalues coincide with the eigenvalues of the 3D hydrogen 
model. 

Now the solutions H„^^^o,o have nonzero total angular momentum while 
zero angular momentum in the 2;-direction, and the logarithm divergence for 
6' = and vr is also present for all ^^^O; indicating that the z-axis plays a 
special role in the classical trajectories analogy. So it is possible to interpret 
that S„£^o,o is related to two-dimensional motions taking place in planes 
containing the z-axis, i.e., to the 2D hydrogen atom. Figure 0] illustrates 
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Figure 3: A boundary surface of \En=i/=ofl{r,9,4>)\'^. 



such interpretation for i = 2,m = 0. The 2D hydrogen atom has also been 
considered in the hterature (see [HUl EH IS21 EBl Ell and references therein), 
but its history is not as controversial as for the ID case. 

Finally, a word about 'En,i,m^o] since they do not belong to the Hilbert 
space, based on the above discussion and proceeding heuristically, it is tempt- 
ing to interpret such solutions as the "contribution" due to the classical tra- 
jectories which come into collision with the nucleus, and the mathematical 
apparatus prudently avoids them explicitly (maybe a mathematical conse- 
quence of the uncertainty principle). 

4 Conclusions 

One is naturally inclined to presume that higher dimensional quantum mod- 
els carry somehow lower dimensional dynamics, and the study of such simpler 
models could mimics important aspects of the original one. Of course, in gen- 
eral the difficulties of performing such dimensional reductions are enormous, 
and usually carried out by "brute force." 

The case of the 3D hydrogen atom discussed in this work has revealed a 
particular and interesting framework: there are experimental evidence for the 
ID hydrogen atom; the 3D hermitian model has just one self-adjoint exten- 
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Figure 4: |^2,o(^)| and |©2,o(^)| (dashed) as function of 6. 



sion, and its 3D eigenvalue equation presents formal solutions q that do 
not belong to the domain of the corresponding Hamiltonian operator; in spite 
of being formal eigenvectors, these solutions live in the underlying Hilbert 
space and present a component in the continuous subspace of the Hamilto- 
nian so that, for an electron in such state, ionization can take place; these 
solutions have formally zero azimuthal angular momentum, with integrable 
probability densities, and are concentrated around the z-axis, indicating their 
ID and 2D character for S„,o,o and 2„£^o,0; respectively. Summing up, such 
solutions are reminiscent of ID and 2D classical trajectories and give a con- 
nection between the hydrogen atom in different dimensions. 

How general is this framework? This is a fascinating open question, whose 
answer could eventually improve the interpretations. 

In addition, notice that there is an attractive relation between the dimen- 
sional interpretations advocated in this work and the mathematical formal- 
ism, which exhausts the possibilities for (pseudo-)eigenvectors. For genuine 
3D motion it presents eigenfunctions in the Hamiltonian operator domain; for 
ID and 2D reminiscent trajectories it presents eigenfunctions in the Hilbert 
space but not in the domain of the operator; and for those colliding trajecto- 
ries (axial divergence) the formal eigenfunctions do not belong to the Hilbert 
space. 
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